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Abstract 

The exact ground state of four electrons in an arbitrary large two leg Hubbard ladder is deduced 
from nine analytic and explicit linear equations. The used procedure is described, and the properties 
^ ', of the ground state are analyzed. The method is based on the construction in r-space of the different 



c/3 



Q 



type of orthogonal basis wave vectors which span the subspace of the Hilbert space containing the 
ground state. In order to do this, we start from the possible microconfigurations of the four 



'^ . particles within the system. These microconfigurations are then rotated, translated and spin- 

O I reversed in order to build up the basis vectors of the problem. A closed system of nine analytic 



linear equations is obtained whose secular equation, by its minimum energy solution, provides the 



^ ■ ground state energy and the ground state wave function of the model. 
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I. INTRODUCTION 



Often, several authors have underhned that in understanding stro ng correlation effects 



presumably new methods of attack and new languages are necessary (jSenthil 



200J). Start- 



ing on this line, in the case of finite value of the concentration of carriers in the studied 
materials, several non-perturbative procedures have been rece ntly worked out , of which 



we mention the fo 
ized in mid-1990s 



lowing ones. 1) Bozonizati on techniques (JTomonagal Il950| ) reactual 



Gulacsi 



( Bowen and Gulacsi 



1997 



Shankai 



1994^ and even extended to finite temperatures 



20011). have been app lied for several sys tems as Coul omb blockade in 



Gulacsi and Anderso: 



Jl998l). 



uantum dots (IGolden and Halperinll2002^ . me soscopic wires (lTorresll2002r). two-band case s 



density waves (jOrignac and Citro 



ferromagn etism ( Honner and Gulacsilll977 ^ 



20031 ). or metal-insulator transitions (JGulacsi and Bedell 



1997 



1998. 19m 



1994^ . The procedure has bee n also used in the s tudy of different model characteristics 



for ex ample, in t he case of t-J (IChen and 



2mM), or Kondo (|McCu11och and Gulacsi 



Wu 



2nm 



McGu11och et a). 



20021). Luther- Emerv (lOrignac and Poiblanc 



2001.20021 ) models. 2) In- 



finite or der canonical trarisformations applied for the Ande rson and multiband Hubbard 



models (IChan and Gulacsi 
Hami 



2001 



Gulacsi and Anderson 



19981 ). and 3) Decompositio n of the 



tonian in positive semidefin ite operators, leading to inte resting new phases ()Gulacsi 



2004a|), al so in three din iensions ([Gulacsi and VoUhardt 



20031 ). or even disordered and in- 



teracting (|Gulacsill2004bl ) cases. 

In the last years, several experimental results have directed the attention to experimen- 
tal situations suggesting interesting technological application possibilities, where the strong 
correlation effects emerge not at finite concentration value of electrons, but for few elec- 
trons confined in a system or device. Such situations related to condensed matter physics. 



are e ncountered on a larg e spectrum of su 



dots (Maksvm et al. 



jfields, as for example i n the case of quantum 



20001). quantum wells (|Kochereshko et al.ll2003l ). mesoscopics ( 



Halfpap 



20011 ). entanglement (jSackett et al.M2000l ). etc. This subject attracted increasing interest 



from the point of view of the theoretic al description in t 



one electron problerns solv ed exa ctlv (ISigrist et al. 



(IKovacs and Gulacsill2001l). three (lA. Amava-Tapia et al. 



1991 



le last decade. Starting from even 
, several cases of interest for two 



2004, 



Dayvdvchev arid Delbourgo 



m 



2004? ) . four (|Zhang and Henlevll2004f ). or few (jPapadopouloi 



2001 



Zhou et al. 



20021) parti- 



cles have been studied. Concretely, in the case of A^p = 4 particles, even if the simulation 



have been started more than a decade ago (|Travnor et al.lll99lL lUesaka et al.lll993h . only 
few valuable results are known in this subject in the condensed matter context, as for exam- 
ple the energy dep endence of the maximal Lyapunov exponent for ID Lenard- Jones system 



Okabe et al.ll2004f) . the b ehaviour of a simplified spinless ferm ion correlated electrons m odel 



Zhang and Henlevll2004[ ). or doped quantum well structures (jKochereshko et al.ll2003h . As 



can be seen, the collected information, at least at the moment, provides only a poor char- 
acterization of the experimental situations of interest. 

In this paper, extending the frame of the non-perturbative methods mentioned above to 
the case of the low density limit, and starting from the aim to provide valuable high quality 
essential information for the Np = 4 case, we provide exact results for this field, presenting 
the exact ground state for four interacting electrons in an arbitrary large two leg Hubbard 
ladder, characterized by periodic boundary conditions. For this to be possible, a direct 
r-space representation is used for the wave functions. On this line first symmetry adapted 
ortho-normalized basis wave vectors are constructed starting from local particle configura- 
tions. Based on these, an explicit and analytic closed system of equations is provided for 
nine type of basis wave vector components, whose secular equation leads to the ground state 
wave function and the ground state energy of the system. 

Deducing the ground state wave function for different microscopic parameters of the 
model, ground state expectation values are calculated for different physical quantities of 
interest, and correlation functions are deduced in order to characterize the ground state 
properties. 

For the case of the low density limit, the procedure presented here, enrols under the 
requirement of new languages. Compairing the deduced ground states to ground states 



obtained in similar conditions for square systems (jKovacs and Gulacsi 



2004f ) at exact level. 



the emerging differences suggest that is highly questionable to approach at a good quality 
level, the two dimensional behaviour from the ladder side. 

The remaining part of the paper is structurated as follows. Section 11. presents the 
Hamiltonian, Section III. describes the used procedure. Section IV. characterizes the de- 
duced ground states, and Section V. containing the summary and conclusions closes the 
presentation. 



II. THE HAMILTONIAN 

The Hamiltonian of the Hubbard ladder is used in a standard form 

H = ti+UU, (1) 

where 

TV 
i= Y. (claCja + CjW), f/ = ^%T%i- (2) 



In these expressions Qo- are canonical Fermi operators which describe electrons on a two 
leg ladder, A^ represents the number of lattice sites in the ladder, {i,j) denotes nearest- 
neighbour sites, t is the hopping matrix element for the electrons, and U is the on-site 
Coulomb repulsion. 

Below we present the exact ground states of the presented model for Np = 4 particles and 
arbitrary ladder lenght, in the case of periodic boundary conditions and u = U/t > 0. The 
deduced ground state is a spin singlet state and from physical point of view describes the 
repulsive Hubbard interaction case (f/ > 0), for t > 0, or the attractive Hubbard interaction 
case (f/ < 0), for t < 0. 

III. THE USED PROCEDURE 

A. The construction of the basis wave vectors 

The procedure we use is as follows. 1) First we number all lattice sites of the ladder 
as shown in Fig.l, taking into account periodic boundary conditions, and considering the 
number of rungs (e.g. A^/2) integer number. 2) We identify the nine possible and quali- 
tatively different microconfigurations in which four particles can be present into the ladder 
(see Fig. 2). At this step we consider that the ladder legs are equivalent, and spin reversed 
microconfigurations are also equivalent. For the clarity of the mathematical notations we 
denote the possible microconfigurations by capital letters A to J, whose lower indices refer 
to the particle positions. In order to clearly distingwish the microconfigurations, we must 
consider in the C,E cases i ^ j, while in the F, J cases j < k (see Fig. 2). We further 
mention that in the process of defining a microconfiguration, at least one of the particles 



is positioned at site number 1. Let us call these particles ,, starter" particles of the micro- 
configuration. For example, for the configuration Ai, the starter particles are two electrons 
with opposite spin; for the microconfiguration Fij^^ the starter particle is one electron with 
spin up, etc. (see Fig. 2). 3) In the following step we define nine different type of orthogonal 
basis wave functions connected to each possible microconfiguration, and denoted with the 
same capital letter introduced in a ket-vector. For example, connected to the microconfigu- 
ration Ai, we have the basis wave vector \Ai); related to the microconfiguration Bi we have 
the basis wave vector |-Bj), etc. 4) Since the lattice sites, legs, and spin orientations are 
equivalent, these properties must be refiected also at the level of the basis vectors. Because 
of this reason, each microconfiguration present in Fig. 2 will be considered as generating 
microconfiguration and denoted by 0\ ■ (n), where n = 1, 2, ..., 9 denotes the microconfigu- 
ration number. For example Of^(l) = Ai,of^(2) = 5^,0^^ (3) = dj, ...^o'l^id) = Jij^k- 
5) Now, connected to each generating microconfiguration 0> (n), we define seven related 



(,, brother") microconfigurations OJ™ (n) (e.g. m = 2,3, ...8) by a) rotating the generating 
microconfiguration m = 1 by 180 degree along the longitudinal symmetry axis of the ladder, 
obtaining the m = 2 brother microconfiguration, b) rotating the generating microconfigura- 
tion m = 1 by 180 degree along the symmetry axis perpendicular to the ladder, obtaining 
the m = 3 brother configuration, c) rotating the m = 3 configuration by 180 degree along 
the longitudinal symmetry axis of the ladder, obtaining the m = 4 brother configuration, 
and finally d) reversing all spin orientations in the m = 1,2,3,4 cases, we obtain the re- 
maining m = 5, 6, 7, 8 brother microconfigurations. In this process it could happen that the 
starter particles situated at site 1 in the generating microconfiguration 771=1, arrive on the 
upper leg for a given brother configuration m > 1. An example of related (,, brother") mi- 
croconfigurations is presented in Fig. 3. for the \Cij) case. 6) At the sixth step, the effective 
construction of the basis vectors follows. A given basis vector \0l ■ {71)), of type n (n being 
considered fixed), connected to the generating microconfiguration OlJ{7i), is constructed 
as follows: a) The different brother microconfigurations OlJ{7i) are all translated along the 
ladder such that the starter particles (see point 1)), arrive on each lattice site of the same 
leg. In this process, the translation must be such effectuated to not modify the interparticle 
positions and relative spin orientations inside the microconfiguration. b) After this step, 
all obtained microconfigurations are added, c) The expression is written in mathematical 
form by representing each microconfiguration by four creation operators acting on the bare 



vacuum with no fermions present. In order to do this, we have to fix the order of creation 
operators for each basis vector type, which has been done as follows. For two doubly occu- 
pied sites, we write the creation operators of the couples next to each other, first the spin 
up, then the spin down contribution as cLcLcLcJi |0). In the case of basis vectors contain- 
ing one doubly occupied site at i we use c||cj|c]|c|,j0). Finally, for basis vectors without 
doubly-occupied sites, the convention cj|c]|c|,|cjj0) is considered, where i < j, k < I must 
hold. 

For example, using the above described steps for the \Cij) basis wave vector, one finds 
the eight related (,, brother") microconfigurations as shown in Fig.3. Furthermore, the math- 
ematical expression of |Cjj) taken for example at i = 2, j = 4, becomes 



1^2,4/ — ( ('^lt'^lj,'^2t'^4i + '^2t'^2i'^3T'^5i + • • 



st At At At 



+ ''^(f +l)T^{f +l)i^{f +2)T^{f +4)i + C(|:+2)T^{f +2)i^(f +3)T^{f +5)i + ■■■'' 
+ (C4|C4|C3|qj^ + C51C51C41C21 + . . .) 

+ *^^(f +4)T^{f +4)i^(f +3)T^(f +l)i + ^(f +5)T^(f +5)i^(f +4)T^(f +2)i + " " ") 
-|-C]^|C^|C4|C2| + C2|C2^C5|C3j^ + ...) + ... j |0) . (3j 

We underline at this step, that because of the fixed conventions presented above, 
somethimes an additional negative sign arises in the process of writting the mathemat- 
ical expressions corresponding to basis vector components shifted from the end to the 
beginning of the ladder in the presence of periodic boundary conditions. For example 
if we shift once more c{|c]Y/2tC2|C3i |0), according to the fixed conventions one obtains 

C2|C;^|C3|C4jO) = — C]^|C2|C3|C4jO). 

B. The deduction of the ground-state wave function 

Our basic observation leading to the solution of the problem is that by applying the 
Hamiltonian H to a basis wave function \OlJ^{ni)) holding a fixed n = rii, only contributions 
of the form \OlJ^J^n)) with n = 1,2, ...,9 can be obtained as a result. Consequently nine 
(e.g. n = 1, 2, .., 9) analytic equations building up a linear system of equations of the form 

H\o^]M) = E E «";/,...|oJ^5',...(^')), (4) 

n'=l i',j',.. 



provide the solution of the problem for arbitrary ladder lenght (e.g. arbitrary A^), where 
a^'J/^ are numerical coefficients. For example, in the case of ra = 1 (e.g. \Ai) = \0l (1))), 
Eq.(jll) becomes 

H\A,) = 2u\A,) - lA,.) - e(z > 2)|Q_i,) - Q{t < iV/4)|Q,+i), (5) 

where Q(K) = 1 if the condition K is satisfied, otherwise Q{K) = 0, and we must have 
1 < 2 < 1 + N/4:. Analogously, forn = 2 (e.g. \Bi) = |o£^ (2))), Eq.© gives 

H\Bi) = 2u\Bi) - Q{t > 1)1 A,.) - e(2 > 1)|A-m) - e(z < Ar/4)|A,m), (6) 

where for the index i one must has 1 < « < 1 + N/A, and similar equations are obtained 
for the remaining n = 3,4,5, ..., 9 index values as well. Consequently, the nine orthogonal 
basis vectors presented in Fig.2. provide nine analytic self-consistent linear equations (e.g. 
Eq.dH)), containing the ground state of the problem. This means that solving the secular 
equation related to Eq.Q we arrive in a subspace of the original Hilbert space which contains 
the ground state. The ground state energy Eg is the minimum possible energy provided 
by the mentioned secular equation, and the ground state wave function is the eigenvector 
corresponding to Eg. The ground state energy and the ground state wave function must be 
numerically obtained from the system of equations Eq.(^. The ground state nature of the 
so obtained solution has been tested by numerical exact diagonalization taken on the full 
Hilbert space at different A^ values. 

The leading terms for two ground state wave functions deduced at a fixed A^ and two 
different u = U/t values are exemplified in the Appendix A. 

The fact that the solution of the ground state of four particles in an arbitrary large two leg 
ladder can be exactly given in Eq.(JH), is related to the observation that the possible n values 
(describing the different type of orthogonal basis vectors entering into the problem) are not 
changing if A^ is increased. Consequently, nine analytical equations will provide the solution 
in Eq. (JU independent of how large the A^ value is. But this does not mean that increasing A^ 
in the course of the numerical treatement of Eq.(j3)), the same number of equations provided 
by Eq.(j3)) are encountered. This is because even if the analytic expression of H\0\ ■ (n)) 
at a fixed n is the same for all A^, the domains covered by the indices i,j, ... in |0j-/ .(^)) 
depend on the A^ value. For example, the Eq.(0) represents the unique analytic equation 
for the \Ai) base vector at arbitrary A^. But during the numerical treatement of Eq.(jSJ, all 



equations for different i values must be considered. Since i < 1 + N/A holds, the unique 
analytic equation (0), will provide an increasing number of different numerical equations by 
increasing A^. But even in this case, the Hilbert space region defined by Eq.(j3]) containing 
the ground state, has an accentuately lower dimension dj.f.fi than the dimension du of the 
full Hilbert space of the problem. For example, at A^ = 16 we have du = 14400, d^ed = 287, 
while at A^ = 32 one has du = 2.410^, but dred = 2141. As seen, at least two orders of 
magnitude reduction in the number of basis wave vectors is encountered in the treatement 
of the problem. 

At the level of principle, Eq.(j3]), based on symmetry properties, delimitates the region 
of the Hilbert space where the ground state is placed. We must however emphasize that 
these symmetry properties, besides the symmetries of the system, depend on the microscopic 
parameters of the Hamiltonian as well. For example, Eq. (j3]) contains the ground state wave 
function only for u > 0. 

IV. THE PROPERTIES OF THE GROUND STATE 

Once the exact ground state wave function is known, the ground state itself can be 
characterized. Starting on this line, in this Section we exemplify the physical properties of 
the deduced ground state at A^ = 28. 

A. Leading terms in the ground state, and ground state expectation values 

We analyze first the u dependence of the leading terms of the ground state wave function 
|\E'g). As seen from Appendix A., at low u the main contributions in l^^g) are obtained 
from relatively closely situated (e.g. almost nearest-neighbour) electrons with opposite spin 
(,,pairs"), the placement of the two pairs being such to maximize the interpair distance. 
Indeed, at u = 3, the leading terms in Appendix A. are the \Dyj), 1-^7,8), IC't.s); l-^e.e) type 
of contributions, which as seen from Fig. 2. (taking into account that the described system is 
a two leg ladder ring with 14 rungs), describe indeed this situation. As long as u is increased, 
the ,, pairs" in the leading terms tend to form double occupancies (as seen in Appendix. A. 
for u = 100), the distance between the pairs remaining considerably high. 

Now we turn to analyze ground state expectation values. The ground state expectation 



value of the kinetic energy term is presented in Fig. 4. As seen, a non-monotonic behaviour is 
obtained. The general tendency for the decrease of the absolute value of Ekm in function of 
u at high u can be understood by the double occupancy formation. The physical reason for 
the presence of the maximum in Ekin for -u < 15 is not yet properly understood. Probably 
the decrease in the increase rate of the double occupancy at a given site above u = 5 causes 
this behaviour (see Fig. 5.). 

The ground state expectation value of the double occupancy per site D = (jii^^ni^i) is 
presented in Fig. 5. The general tendency present in this figure is that the increasing u value 
leads to the increase of D. This can be understood if we remember that for t < the 
analyzed situation describes the attractive on-site interaction case, hence by increasing u we 
increase the number of double occupancies. 

B. Pair correlation functions 

Pair correlations are analyzed via the density-density correlation function 

1 ^ 
C„(r) = —Y^{{nihi^.,) - {ni){ni+r)), (7) 

^^ i=i 

and spin-spin correlation function defined by 

-, N 1 -^ 

4 = 1 j = l 

where hi = h^ + hii, hi^„ = cl^Cif^, Sf = (hi^^ — hi^i)/2, r measures the inter-site distance in 
lattice constant units, and (...) has the meaning of the ground state expectation value. 

The behaviour of the spin-spin correlation function is presented in Fig. 6. As seen, com- 
paired to the non-interacting case, the spin-spin correlation decreases if u is increased. Con- 
cerning the distance dependence of Cs'{r) at a fixed u, we see that even the short-range 
correlations are strongly (presumably exponentially) decreasing, and the spin correlation 
lenght covers practically only nearest-neighbour sites. 

The density-density correlations are exemplified in Fig. 7. The presented curves have a 
specific structure which can be understood based on the analyzes of the leading terms of the 
ground state wave function presented at the beginning of this Section. As observed there, 
in the ground state, two ,, pairs" tend to be situated at highest possible distance each from 
other, providing the behaviour presented in Fig. 7. 



V. SUMMARY AND CONCLUSIONS 

We have deduced the exact ground state for four electrons in an arbitrary large two-leg 
Hubbard ladder and have analyzed the physical properties of the ground state in function 
of microscopic parameters of the model. The procedure is based on the construction in 
r-space of nine different type of orthogonal basis vectors which span the subspace of the 
Hilbert space containing the ground state. In order to do this, we start from the possible 
microconfigurations of the four particles within the system. This microconfigurations are 
then rotated, translated and spin-reversed in order to build up the basis vectors of the 
problem. A closed system of linear equations is obtained whose secular equation, by its 
minimum energy solution, provides the ground state energy and the ground state wave 
function of the model. The dimensionality of the subspace containing the ground state is 
substantially less than the dimension of the full Hilbert space. The deduced ground state 
wave functions have been used for the calculation of ground state expectation values and 
correlation functions in the process of the characterization of ground state properties. 
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APPENDIX A: EXPLICIT GROUND STATE WAVE FUNCTIONS. 

We present below the leading terms of explicit ground state wave functions deduced for 
A^ = 28, at -u = 3 and u = 100. The ground state |\l'g) is normalized to unity, and contains 
ortho-normalized basis vectors. 

At -u = 3 we obtain for the ground state wave function 

\^,iu = 3.0)) = 

0.l776U\Drj) + 0.1700071^7,8) + 0.1 70006 1^7,8) + 0.158781^6,6) 
+0.158091|C6,7) + 0.1580891^6,7) + 0.138593|D8,7) + 0.138593|L'7,8) 
+0.135254|C5^6) + 0.1352361^5,6) + 0.1291991^5,5) + 0.1288691^6,7) 

10 



+0.128868|L>7,6) + 0.110223|L>5,6) + 0.110212|L)6,5) + 0.103555|C4^5) 
+0.1035021^4,5) + 0.100914|G7,8,i) - 0.100914|G'8,2,8) + 0.09744121^6,8) 
+0.0974383|C6,8) - 0.09383711^7,2,7) + 0.09383651^6,7,1) + 0.093035|D7,9) 
+0.0917791 |L>4,4) + 0.08974591^8,6) + 0.0897451^6,8) + 0.0870033|C5,7) 
+0.08699741^5,7) + 0.0843159|L'4,5) + 0.0842604|L)5,4) - 0.0822647| J7,i,8) 
+0.082264|/J7,7,i4) - 0.08026861^6,2,6) + 0.08026251^5,6,1) + 0.0801191^5,7) 
+0.08011231^7,5) + 0.07743511^7,8,14) - 0.07718431^8,2,9) + 0.0764892|i/6,6,i4) 
-0.07648891 J6,i,7) + 0.07469791^6,7,14) - 0.07469781^8,2,7) + 0.0744824|G7,8,2) 
-0.0744821 |G7,2,8) + 0.07049931^4,6) + 0.0704551^4,6) + 0.06668691^5,6,14) 
-0.06668611^7,2,6) - 0.0666109|G6,2,7) + 0.06660521^6,7,2) + 0.0663597|C3,4) 
+0.06631151^3,4) + 0.0654117|iJ5,5,i4) - 0.0654063| J5,i,6) + 0.0650782|L)8,8) 
+0.064882|D4,6) + 0.06483051^6,4) - 0.06142561^5,2,5) + 0.06139381^4,5,1) 
+0.05993051^7,9,1) + 0.05946961^6,9) + 0.0594586|C6,9) + 0.05781231^6,8,1) 
-0.05781 18 1 G8,3,8) + 0.05747751^7,10) + 0.05747451^6,9) - 0.0555399|if6,i3,7) 
+... 

while for u = 100 one has 

l^giU = 100)) = 

0.384661^7) + 0.3846111^7) + 0.3547281^6) + 0.3545581^6) 
+0.3066861^5) + 0.3060651^5) + 0.243462|54) + 0.2411791^4) 
+0.1699811^3) + 0.1615761^3) + 0.139591^8) + 0.1395811^8) 
+0.0961638|S2) + 0.06524681^2) + 0.06235881^7,8) + 0.0623528|C7,8) 
+0.06156561^7,7) + 0.05915131^6,7) + 0.0591339|C6,7) + 0.05676491^6,6) 
+0.05291331^5,6) + 0.0528502|C5,6) + 0.0490391 |D5,5) + 0.0440121^4,5) 
+0.0437798|C4,5) + 0.03878631^4,4) + 0.03307531^3,4) + 0.0322208|C3,4) 
+0.026534811)3,3) + 0.0223424|D8,8) + 0.02129071^2,3) + 0.018145|C2,3) 
+0.0171881|5i) + 0.01290861^2,2) + 0.011569|Ei,2) + 0.00249322|D8,7) 
+0.00249322|D7,8) + 0.002364741^6,7) + 0.00236474|L)7,6) + 0.00211441|L)5,6) 
+0.00211441|D6,5) + 0.00175512|L)4,5) + 0.00175512|D5,4) + 0.00130539|L'3,4) 

11 



(Al) 



+0.00130539|L>4,3) - 0.0012476|G8,2,8) + 0.0012476|G7,8,i) + 0.001215591^6,8) 
+0.00121536|C6,8) + 0.00118331|G6,7,i) - 0.00118331|G'7,2,7) + 0.0011211|E5j) 
+0.00112029|(:75 7) - 0.001058051^6,2,6) + 0.001058051^5,6,1) + 0.0009696431^4,6) 
+0.000966692|C4,6) - 0.000878259|G'5,2,5) + 0.0008782591^4,5,1) + 0.000788073|L>2,3) 
+0.000788031^3,2) + 0.0007711791^3,5) + 0.000760318|C3 5) - 0.000653214|G4,2,4) 
+.... (A2) 
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FIG. 1: The numbering of the lattice sites for the two leg ladder taken with periodic boundary 
conditions. N is considered even. 
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FIG. 2: The different possible types of basis vectors. We note that for the cases C,E i ^ j, 
while for F, J j < k is considered, respectively. In the cases F, G, H, J, the double occupancy is 
forbidden. 
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FIG. 4: The dependence of the kinetic energy on u. 
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FIG. 5: The dependence of the double occupancy per site D = |--Ty^| on 
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FIG. 6: The density-density correlation function for n = (dots, dot-dashed hne), u = 10 (squares, 
long dashed line), u = 30 (diamonds, short dashed line), u = 100 (stars, continuous line), r is the 
distance in lattice constant units 
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FIG. 7: The 5^-5*^ correlation function for u = (dots, dot-dashed hne), u = 10 (squares, long 
dashed line), u = 30 (diamonds, short dashed line), u = 100 (stars, continuous line), r is the 
distance in lattice constant units 
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